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Abstract 

A spin-isospin dependent Three-Dimensional approach based on momentum vectors for formu- 
lation of the three-nucleon bound state is presented in this paper. The three-nucleon Faddeev 
equations with two-nucleon interactions are formulated as a function of vector Jacobi momenta, 
specifically the magnitudes of the momenta and the angle between them with the inclusion of the 
spin-isospin quantum numbers, without employing a partial wave decomposition. As an applica- 
tion the spin-isospin dependent Faddeev integral equations are solved with Bonn-B potential. Our 
result for the Triton binding energy with the value of —8.152 MeV is in good agreement with the 
achievements of the other partial wave based methods. 
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I. INTRODUCTION 



During the past years, several methods have been developed to solve the nonrelativistic 
Schrodinger equation accurately for few-nucleon bound states, by using realistic nuclear 
potentials. These methods are the CRCGV fl, the SV Q, the HH Q, the GFMC £|, the 
NCSM Q, the EIHH |g] and the Faddeev. These calculational approaches are mostly based 
on a partial wave (PW) decomposition. Stochastic and Monte Carlo methods, however, are 
performed directly using the position vectors in the configuration space. One of the most 
viable approaches appears to be the Faddeev method. 

The calculations based on the Faddeev approach are performed after a PW expansion with 



phenomeno 
space [T~ 



ogical potentials either in the momentum space [7|- 1^] or in the configuration 



18l | . Recent bound state calculations with the Faddeev approach have been done 



with the chiral potentials in the momentum space 19|-[22j. Experience in three- nucleon 
calculations shows that the standard treatment based on a PW decomposition is quite 
successful but also rather complex, since each building block related to involved operators 
requires extended algebra. The Faddeev calculations based on a PW decomposition, which 
includes the spin-isospin degrees of freedom, after truncation leads to a set of a finite number 
of coupled equations in two variables for the amplitudes and one needs a large number of 
partial waves to get converged results. In view of this large number of interfering terms it 
appears natural to give up such an expansion and work directly with the vector variables. 
On this basis three- and four-body bound states have recently been studied in a Three- 
Dimensional (3D) approach where the spin-isospin degrees of freedom have been neglected 



in the first attempt 



23]- [27|. In the case of three-body bound state the Faddeev equations 



have been formulated for three identical bosons as a function of vector Jacobi momenta, 
with the specific stress upon the magnitudes of the momenta and the angle between them. 
Adding the spin-isospin to the 3D formalism is a major additional task, which will increase 
more degrees of freedom into the states and therefore will lead to a strictly finite number of 
coupled equations [28| . In this paper we have attempted to implement this task by including 
the spin-isospin degrees of freedom in the 3N bound state formalism. To this end we have 
formulated the Faddeev equations for the 3N bound state with the advantage of using the 
realistic NN forces. The presented 3D formalism in this paper in comparison with the 
traditional PW formalism avoids the highly involved angular momentum algebra occurring 
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for the permutation operators. According to the spin-isospin states that have been taken 
into account, we have obtained the eight, twelve, sixteen and twenty four coupled equations 
for a description of the 3N bound state, i.e. 3 H and 3 He. In this way, we solve the Faddeev 
integral equations for calculation of the Triton binding energy with Bonn-B potential. The 
input to our calculations is the two-body t-matrix which has been calculated in an approach 
based on a Helicity representation and depends on the magnitudes of the initial and final 



29]. 



momenta and the angle between them 

This manuscript is organized as follows. In section [III we present the formalism. Meaning 
that we have derived the Faddeev equations and the 3N wave function in a realistic 3D scheme 
both as a function of Jacobi momenta vectors and the spin-isospin quantum numbers. Also 
the novel 3D representation of the Faddeev equations is contrasted with the corresponding 
traditional PW representation. In section II III we present our results for the Triton binding 
energy and compare them with the results obtained from the PW calculations. In order 
to test our calculations the calculated expectation values of the Hamiltonian operator are 
compared to the obtained eigenvalue energies. Finally in section [IV] a summary and an 
outlook will be presented. 



II. FORMULATION FOR 3N BOUND STATE IN A 3D FADDEEV SCHEME 
A. The Faddeev Equations 

The bound state of three pairwise-interacting nucleons is described by the Faddeev equa- 
tion jllj: 

|V> = \ik 2>3 ) = G tP\i(>), (1) 

where Go is the free 3N propagator, t denotes the NN transition matrix determined by a 
two-body Lippman-Schwinger equation and P = P12-P23 + P13P23 is the sum of a cyclic and 
anti-cyclic permutations of the three nucleons. The total 3N wave function \*f?) is composed 
of the three Faddeev components as: 

|¥) = (1 + P)M. (2) 

The antisymmetry property of \ifj) under exchange of the interacting particles 1 and 2 
guarantees that |^) is totally antisymmetric. In order to solve Eq. ([1]) in the momentum 
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space we introduce the 3N basis states in a 3D formalism as, (see Fig. [[]): 



|pqa) = \pqa s a T ), (3) 

the basis states involve two standard Jacobi momenta p and q [ll], and | a) is the spin- 
isospin parts of the basis states, where the spin part is defined as: 



a s ) = |((si s 2 )s 12 s 3 )SMs) = \(s 12 -)SM S ), 



(4) 



and the isospin part | «t) is similar to the spin part. As indicated in Fig. [T] the angular 
dependence explicitly appears in the Jacobi vector variables, whereas in a standard PW 
approach the angular dependence leads to two orbital angular momentum quantum num- 
bers, i.e. lu and Z3 llj. It indicates that in the present 3D formalism there is not any 



coupling between the orbital angular momenta and the corresponding spin quantum num- 
bers. Therefore we couple the spin quantum numbers su and S3 to the total spin S and its 
third component Ms as: | (s 12 s^S Ms ). For the isospin quantum numbers similar coupling 
scheme leads to the total isospin T, Mt as | (£12 t 3 )T Mt). 



3D 



PW 





j3 y 
1 *i 

|p q (s 12 1/2)S M s (t 12 1/2)T M T > |p q ((l 12 s 12 )j 12 (l 3 l/2)j 3 ) J Mj(t 12 1/2)T M T > 



FIG. 1: Definition of the 3N basis states in the 3D approach in comparison with the corresponding 
basis states in the PW approach. 



In order to evaluate the transition and the permutation operators we need the free 3N 
basis states | p q 7 ), where 



i) = \1s1t), lis) = \m sl m S2 m SA ). 



(5) 



The quantities m Si (i = 1,2,3) are the third components of the spins of the three nucleons. 
The isospin part of the basis states \jt) is similar to the spin part. To achieve this aim 
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when changing the 3N basis states | a ) to the free 3N basis states | 7 ) we need to calculate 
the following Clebsch-Gordan coefficients (see appendix |A|) : 

(7l«) = 9^ = 9^9^ = (m Sl m S2 m S3 \(s 12 -)S M s ) (m tl m t2 m t3 \(t 12 -)TM T ). (6) 

The introduced basis states are complete and normalized as: 

Y,jd?p J ^|pq£)(pq£| = l, (pq?|p , q / r) = 5 3 (p-p / )5 3 (q-q , )^ f) (7) 

where £ indicates a and 7 quantum number sets. Now we can represent the Eq. (CQ) with 
respect to the basis states which have been already introduced in Eq. ([3]): 

(pqa|V>) = Yl J d3p> / d 3 q' (p qa\G tP\ p'c^ a') (p'c^ a' \ip). (8) 

a' 

For evaluating the Eq. ([8]), we need to evaluate the matrix elements of 
( p qa \GotP\ p' q' oi ), towards this aim, it is convenient to insert the free 3N complete- 
ness relations as: 

(pqa\G tP\p q a) = ^ E ( a I 7 )( P q7 \G tP\ p q' 7 )( 7 | a' ) 

7 7' 

= X] &«7 07'<*'(P < 17|G ! O*-P| P'q'Y)- (9) 

7:7' 

For evaluating the matrix elements ( p q7 \ GotP\ p' q' 7' ) we should insert again a free 3N 
completeness relation between the between the two-nucleon t-matrix and the permutation 
operators as: 

(PC^IGO^IPVY) = E _£_3£ E / ^P" J ^l" 

m 4m, 7" 

x (pqtrp q 7 )(p q 7 \p\ p q 7 ), (io) 



where the matrix elements of the two-body t-matrix and the permutation operator P are 
evaluated separately as: 

(pq7|t|p"q ,/ 7 ,/ ) = 5 3 (q - q") 5 Ws3 m» 3 ( P™^™^ m tl m t2 |i(e)|p" m*m* m'^m'Q, 

(11) 

(p-q'vi^ipvy) 

13 1 

= ^(p" + 2? + 4^') ^fa" - P' + ^<l<2 5 <' 2 < 3 ^i'sWii 5 K>t 2 <W 2 <, <W 3 K a 

13., 1 

+<5 (P + 2? ~~ 4^ ) ^ + P + 2^ ) ^m^m^ dm'J 2 m' sl $m'J 3 m' S2 ^m'^m'^ ^m^m^ ^m^m^, 



(12) 



where the two-body subsystem energy in the NN t-matrix is e = E — I 2 -. 

In order to evaluate the matrix elements of the permutation operator P we have used 
the relation between the Jacobi momenta in the different 3N systems (312), (231) and (123). 
Inserting Eqs. ([II]) and (H2D into Eq. (JTOD leads to: 

{VQl\GotP\ p'q'V) 



J? P 2 3g 2 
m Am 



x <j ^ 3 (q - P' + ^q') $m S3 m> sl $mt 3 m' H (pm Sl m S2 m tl m t2 \t(e)\^-q - q'm' Sa m' Sa m' t2 m' t3 ) 

+5 3 (q + p' + -q') 5 m33 m' S2 ^m t3 m; 2 (pm Sl m S2 m tl m <2 |*(e)|~q + q m s .rd Si m' ta m' tl ) 

(13) 

Inserting Eq. ([TBI into Eq. Q and consequently inserting into Eq. §E§ and integrating 
over p' variable yields: 



(pq«M = ^_ g _ gg2 9^ gy a > J d 3 q' 

m Arm 'v 'v' /-v' 



m 4m, 7, 7', a' 

x (pm sl m S2 m tl m t2 |t(e)|^q - q' m' S2 m' S3 m' t2 m' t3 ) 5 ms3<i Smt 3 m' H (q + ~q' q' a'l^) 

+ (pm Sl m S2 m tl m t2 |t(e)|^q + q , m / S3 m / Si m; 3 m; i ) <W S3 < 2 <W 3 mJ 2 (-q - ^q' q'a'|<0) 

(14) 

Applying the permutation operator P\ 2 action on the Faddeev component, the space and 
also the spin-isospin parts of the basis states, results in: 



PM) = 


-\1>), 




ftalpq) = 


1 -pq)> 




Pn\a) = 






Pn\l) = 


m S2 m Sl m S:j m t2 m tl m t3 ) , 


(15) 



and consequently the following relations would be concluded: 

(pqa|^) = -(-r i2+tl2 (-pqa|V>, 
(pm Sl m S2 m^mtal^lp'm^m^ m' tl m' t2 ) = (pm Sl m S2 m tl m t2 \t(e) P 12 \-p' m' S2 m' sl m' t2 m' h ). 

(16) 



Therefore, we can rewrite Eq. tHM as: 
ipq«l^) = E _ g _ 3^ Yl 9<x 7 9ya> J d 3 q 



m 4m 7, 7', a' 

x { {vm Sl m S2 m tl m t2 \t(e)\^^ - c(m' S2 m' ss m^ 

x(q+ ^q' q'a'l^) 
+ (pm Sl m S2 m tl m i2 |t(e) Pi 2 |^q - q'm^m^ m^m^) 5 ms3<2 <W 3 m' t2 



;_)s 12 +* 12 )( q+ _ q 'q' a 



p 2 3g 2 ^ ^ 9ay 9~/'a' ^rn ss m' ai $mt :j m' ti J d q 

n-i /I -wi f\/ f\l' r~v? 



E 

m Am 7,7 ,a 



x (pm sl m S2 m tl m t2 |i(e)(l - Pi 2 )|— -q- c( m' S2 m' S3 m' t2 m' t3 ) (q+^q'q'a" 



(17) 



The final derivation of Eq. ffl7]) is made by the exchange of labels m' , m' tl to m' S2 , m' t2 
and reverse of it in the second term as well as the following relation; 

Of* = R S ' 12+t ' 12 ( < < < I ( s 'i2 \)S' M' s ) ( m' t2 m' tl m' ta \ (t' 12 ~)T' M' T ) . (18) 
By introducing the physical representation of the two-body t-matrix follows (see appendix 



a (pm Sl m S2 m tl m t2 \t(e)\ p m' si m' S2 m' tl m' t2 ) a 

= (pm Sl m S2 m tl m t2 \t(e)(l - P l2 )\p m Si m S2 m' ti m t2 }, (19) 

the three-dimensional Faddeev integral equations can be obtained as: 

(pqa \lp) = — J 2 7^2 ^ gar f g ~l' a ' ^ m *3<l ^3™^ 

m 4m 7, 7', a' 

f — 1 1 

x / d 3 q a (pm Sl m S2 m tl m t2 \t(e)\—q_-^m' S2 m' Ss m[ 2 m' t3 ) a (q + -q' q'a'l^)- (20) 

The Faddeev component (pqa|i/>) is given as a function of Jacobi momenta vectors, 
p and q, and also quantum number sets, Ot . clS cl solution of the spatial three-dimensional 
integral equations, Eq. (120]) . In order to solve this equation directly and without employing 
the PW projections, we have to define a coordinate system. It is convenient to choose the 
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spin polarization direction parallel to the z-axis and express the momentum vectors in this 
coordinate system. By these considerations we can rewrite Eq. fl20l) as: 

1 



^ a (PQX P q) = 3 ^2 I dq'q' 2 I dx q > I d(f q 



E- 

m 4m 



dq' q' 2 / dx q i I 

I! J-l JO 



p 2 3g 2 

m 4m 

; e) ijj a (it 4 Xnq'), (21) 

a 1 

where 



/ m i m / 



~ \ V — * <• _ S2 s 3 t 2 f 3 

T aa '{p, 7T, a^pff; e) = y ^ 9ay 9ya' Um sz m' Sl ^mt :j m' ti ta m ai m S2 rnt Y mt 2 \P-i 7I ") ^pfrj e ) ; (22) 



7.7' 



%>g — ^p^g + -w 1 — ^t,-\/1 — sin((/)„ — 



= x p x ql + J I - J I - x 2 , sin(0 p 



Xqq> = X q X q t + J 1 - X 2 \ / 1 - XjL Sm(<^ 



^ = xhv 2 + 4 2 + qq'xqq', 



Xpn 



1 

■^qXpq ~\- q Xpqi 
7V 



1 



71 = \l q 2 + ^4 2 + qq'xqq', 

71 

In a standard PW approach, Eq. ( |2TI) is replaced by a set of an infinite number of coupled 
two-dimensional integral equations for the amplitudes with the kernels containing relatively 
complicated geometrical expressions: 



1 



f (p ?) = - — -2 — ^2 / dq'q / dx q > 



+1 



V ttt^ G aa >(q,q,x q >) p , (24) 



7T'12 7T'i2 



where, as is shown in Fig. [TJ the spin-space as well as the isospin parts of the basis states 
in the PW decomposition are \a) = |((Zi2Si 2 )ji2 [hSci)h)JMj (t 12 t 3 )TM T ). G aa >(q,q',Xq') 
is composed of Legendre functions, powers of q and q' an d purely complicated geometrical 
quantities like Clebsch-Gordan coefficients and Qj symbols. The comparison of Eqs. ( 12T]) and 
(1241) shows that new 3D formalism avoids the highly involved angular momentum algebra 
occurring for the per mutations and additionally it will be more efficient especially for the 



three-body forces 



48) 



B. The 3N Wave Function 

The representation of the total wave function, Eq. (j2j), with respect to the basis states 
which have been introduced in Eq. ([3]), reads as follows: 

(pqa|*) = (pqa|(l + P)|V>) 

= (pq a\i/}) + (pq a|P 12 P 23 |<0) + ( p q a |P 13 P 23 |^>, (25) 

where the first Faddeev component 

(pqa\ip) = 3 (pqa|^) = f (p,q) =f (PP M ), (26) 

is given explicitly as a three-dimensional integral equation, Eq. (1211) . Here the subscript 
3 of the bra basis states stands for the three-body subsystem (12,3), which as matter of 
convenience, is called subsystem 3. For the second and third components we need to evaluate 
the action of the cyclic and the anti-cyclic permutation operators P12P23 and Pi 3 P23 on the 
first component as: 

(pq a IP12P23IV') 



(pq a IP13P23IV') 



the space as well as the spin-isospin parts of the coordinate transformations 3(1)1 and 3(1)2 
can be evaluated as: 

3 (pqa|p'q / a / ) 1 = 3 ( p q | p' q' ) 1 3 ( a | a' ) 1 

= S 3 {p' + ^p+-q)6 3 (q[-p + -q) 

X $MsM' s $SS' &M T M' T $TT' C* S {(X S , s' 23 ) C^(a T , t' 23 ) , ( 28 ) 



3 (pq a (P^Pal^) 

J d 3 p' J d 3 q' 3 (pqa\P 12 P 23 \p' ^ a' ) 3 3 (p' ^ a' \ 

a' 

Y, f d 3 p' f dV 3 (pq«|p / q / « , )i3(p , q / « , |^), 



3 (pq a |Pi 3 P 23 |^) 

J d 3 P ' f d 3 q' 3 (p q a\P 13 P 23 \p'^a') 33 (p'^a'\ij) 

a' 

Y J d^V f d 3 q 3 (pqa\ p'q / a / ) 2 3(p'q'a / W), (27) 
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3 (pqa| p / q / a) 2 = 3 ( P Q I p' q' >2 3 ( « I «' )2 

= 5 3 (p' + ^p-Jjq)W + p + ^q) 

X <VsMj, <W &M T M' T $TT' Cg*(a s , S31) Qr*( a T, £31), (29) 

where the spin coefficients C£ and C£* are given as: 

QK,4 3 ) = (-)^ + ^ 1+S2+S 3 I *i *» ^ 

I s 3 S s' 23 

CF(«5,4) = (-) s 3 1 +2 S2+S3+ s 1 J S l ^ ^12 . (3Q) 



s 3 S si 



31 



and the isospin coefficients and C?* are similar to the corresponding spin coefficients. 
By these considerations we obtain the second and third Faddeev components as: 



(pq a |P 12 P 23 |^) = J2 C5(«5,4)^TK,4)^(-^P-|q,P-^q) 

s' t' 

= E c s( a s, s'23) Cr(oiT, ^23) ^'(^i ^2 x Wl7T2 ), 

a 23> 23 

(pqa|P 13 P 2 3|V) = E C 5*K,4i)CrK,4)^(-^P+^q,-p-^q) 

*31' , '31 

= E ^(a s y 31 )C* r *(a r) 4)f "(n! n 2 xn.nj, (31) 



s 31' t 31 



where 



«*> = 1(4 ^M 5 (4 i)TM T >, 



'1 „ 9 „ 3 

^ = \l-^ 2 ^YQ q +l pqXpq ' 



K-2 = \/p 2 + ^q 2 -pqx pq , 



= -^—{-\p 2 + \q 2 -\p ( l X pq), (32) 
7Ti7T 2 Z O Z 
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1 ** \ 

a ) = 


1 ( s si ^)SM S (t u -)TM T ) 




I 1 2_L 9 2 3 


III = 




n 2 = 


\jv 2 + ^Q 2 +P(ix pq , 




/ 2 ^ 2 I > 


^UiU 2 = 


(-» PQXva 

nin 2 v 2 F 8 y M ' 



(33) 



C. Comparison of Coupled Faddeev Equations in both 3D and PW Schemes 

In this section we discuss the number of coupled equations in both 3D and PW approaches. 
In a standard PW approach the infinite set of coupled integral equations, given in Eq. (1241) , 
is truncated in the actual calculations at sufficiently high values of the angular momentum 
quantum numbers. If one assumes that the NN t-matrix acts only in very few partial waves 
then the number of the coupled equations are correspondingly small. As shown in Table [J 



if NN t-matrix acts up to j^ ax = 1,2,3,4 and 5, then the number o 
26, 34 and 42. This is while the total isospin is restricted to T = \ 



' channels will be 5, 18, 
321. 



TABLE I: The number of PW channels which compose the Triton wave function when the NN 
i-matrix acts up to different total two-nucleon angular momenta jyf x . Total isospin is restricted 
to T = |. The number of channels for j^ ax = 1, namely iV a = 5, is related to only positive parity 
states. 



■max 
J12 


1 


2 


3 


4 


5 


N a 


5 


18 


26 


34 


42 



In Table [TT] we list all the spin-isospin states which compose the 3N, i.e. 3 H and 3 He, 
wave function and consequently in Tables II III and IIVI we present the number of spin-isospin 
states for the 3N bound states as well as the number of coupled Faddeev equations in 
realistic 3D formalism presented in this paper. It is clear that = 4r refers to 3 He and 
Mt = ~y refers to 3 H. Since the angular momentum quantum numbers, i.e. lu,h, do 
not appear explicitly in our formalism, therefore the number of coupled equations which are 
fixed according to the spin-isospin states are strongly reduced. This is an indication that the 
present formalism automatically considers all partial waves without any truncation on the 
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TABLE II: Quantum numbers of the spin-isospin states which compose 3 H or 3 He wave function. 



channel 


(S12 \)S M S 


(ii9 


\)TM T 


fs 


- T) 


1 


m hi 

(° 2)2 


+1 
2 


fo V 

V u 2- 


1 
2 


+1 /-1 
2/2 


^2 


- h) 

2! 


2 


(° 2)2 


-1 
2 


l u 2- 


1 
2 


+1 /-1 

2/2 


V2 


2> 


3 


1 1 1 \ 1 
(! 3)5 


ir 


fo V 

2 / 


1 
2 


+1 / — 1 
"2 / "2" 


V2 


2> 


A 


I 1 2)2 


-1 
2 


fo i N 


1 
2 


+1 /-1 
2 / 2 


1.2 


- ±1 
2^ 


5 


(° 2)2 


+1 
2 


(1 ^ 

V 2' 


1 
2 


+1 /-1 
2 / 2 


V2 


2/ 


6 


(° 2)2 


-1 
2 


(1 ^ 

V x 2- 


1 
2 


+1 /-1 
2/2 


V2 


2> 


7 


/1 1 \ 1 


+1 
T" 


(1 l" 

V 2' 


1 

2 


+1 / -i 
"2"/ "2" 


(h 

\2 


2> 


8 

O 


I 1 2)2 


-1 
2 


fl i % 
V 1 2- 


1 
2 


+1 /-1 

2 / 2 


fl 


2^ 


9 


(° 2)2 


+1 
2 


(1 V 


3 
2 


+1 /-1 

2 / 2 


V2 


" 1) 

2/ 


10 


/ n hi 
1° 2)2 


-1 
2 


(1 V 

V 2' 


3 
2 


+1 /-1 

2 / 2 


V2 


2/ 


11 


I 1 2)2 


+1 
2 


(1 V 

V x 2^ 


3 
2 


+1 /-1 
2 / 2 


V2 


2> 


12 


I 1 2)2 


-1 
2 


(1 ^ 

V 2' 


3 
2 


+1 /-1 
2 / 2 


(h 

V2 


-|) 

2/ 


13 


I 1 2)2 


+3 
2 


fo i x 


1 
2 


+1 /-1 
2/2 


f 5 


2> 


14 


I 1 2)2 


+1 
2 


fo i x 

V u 2< 


1 
2 


+1 /-1 
2 / 2 


V2 


- l ) 

2> 


15 


(! 5)1 


_i 
T" 


fo V 

V u 2 / 


1 

2 


+1 / — 1 
"2/ "2" 


fi 
^2 


2> 




I 1 2)2 


-3 
2 


fo i % 

l U 2- 


1 
2 


+1 /-1 
2 / 2 


f 3 - 

V2 


- ±1 
2> 


17 


I 1 2)2 


+3 
2 


fl V 

V 2' 


1 
2 


+1 /-1 
2 / 2 


(1 
V2 


-h) 

2> 


18 


n 1 \ 3 
I 1 2)2 


+1 
2 


fl V 

V 2' 


1 
2 


+1 /-1 

2 / 2 


V2 


- h) 

2> 


19 


u 2 J 2 


-1 
2 


(1 V 

V x 2' 


1 
2 


2 / 2 


fl 
V2 


- h) 

21 


20 


a m 


-3 
2 


v 2 ' 


1 
2 


±I/=I 

2 ' 2 


(1 
v 2 


-5) 


21 


a £)§ 


+3 
2 


a 


3 
2 


2 / 2 


V2 


2^ 


22 


a £)§ 


+1 
2 


a I 


3 
2 


2 / 2 


V2 


2^ 


23 




-1 
2 


a 2: 


3 
2 


2 / 2 


^2 


2^ 


24 




-3 
2 


a 2: 


3 
2 


2 / 2 


( 5 
V2 


2) 



space part. Considering the spin-isospin degrees of freedom for both 3 H and 3 He states yields 
the same number of coupled equations and it leads to 8, 12, 16 and 24 coupled equations 
for different combinations of the total spin-isospin states S — T: (| — |), (| — |*), (|* — |) 
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TABLE III: The number of spin-isospin states for 3N bound sates, i.e. 3 H and 3 He, in a realistic 
3D formalism. N$ and Nt are the number of spin and isospin states correspondingly. 



(812 h) S M S 


s = i 


s = l 


c 13 
° 2> 2 


(*12 i)T M T 


- 1 2 


rp 3 

J — 2 


T 13 
- 1 2' 2 


(Ol)IM 


2 





2+0 


2)2 ^2~/ "IT 


1 





1+0 


fl 1)1 ±i 

K 1 2)2 2 


2 





2+0 


fl 1)1 ±1/^1 
V 1 - 2/2 2 / 2 


1 





1+0 


(-i 1\3 ±1 ±3 
V 1 2)2 2 2 





4 


0+4 


fl 1)3 ±1/^1 
v x 2/2 2 / 2 





1 


0+1 


N s 


4 


4 


8 




2 


1 


3 



TABLE IV: The number of coupled Faddeev equations for the 3N bound state, i.e. 3 H and 3 He, 
in a realistic 3D formalism according to the spin-isospin states (S — T). N = N$ x Nt is the total 
number of coupled Faddeev equations. The star superscript indicates all the spin or isospin states 
that one can take into account up to a specific value. 



(S-T) 


(l 1\ (l 3*\ /3* 1 
V2 2) \2 2 ) \2 2 


) (I*"!*) 


N s 


4 4 8 


8 


N T 


2 3 2 


3 


N 


8 12 16 


24 



and (|* — I*) respectively. The star superscript indicates all the spin or isospin states that 
we have taken into account up to a specific value. It is clear that in the 3D formalism, e.g. 
for a fully charge dependent calculation, there is only 24 coupled equations, whereas in the 
PW approach after truncation of the Hilbert space to T — | there is 42 coupled equations. 
Therefore our 3D formalism leads to a small number of coupled equations in comparison 
with the very large number of coupled equations in the truncated PW formalism. However, 
it should be mentioned that our formulation leads to coupled equations in three variables 
for the amplitudes, whereas the PW formulation after truncation leads to a finite number 
of coupled equations in two variables for the amplitudes. So the 3D formulation leads to a 
lesser number of coupled integral equations in three dimensions and the PW formulations 
leads to more coupled integral equations in two dimensions. Thus, the price for the smaller 
number of equations is the higher dimensionality of the integral equations. In other words, 
algebraic simplification is achieved by a more involved numerical scheme. 



13 



III. NUMERICAL RESULTS FOR 3 H 



A. Triton Binding Energy 

In order to be able to test our realistic 3D formalism for the 3N bound state we solve 
the three-dimensional Faddeev integral equations, Eq. (|2ip . We calculate the Triton binding 
energy by solving eight coupled Faddeev equations for (| — ~) spin-isospin states and compare 
our results with the other PW results. In this res pec t, we use Bonn one-boson-exchange 
(OBE) potential in the parametrization o: 



Bonn-B 32j and in an operator form which can 
be incorporated in the 3D formalism [29|. In the numerical treatment, the dependence 
of Faddeev components to the continuous momentum and the angle variables, should be 
replaced by a dependence on certain discrete values. For this purpose we use the Gaussian 
quadrature grid points. 

TABLE V: The calculated binding energy Et of the three-dimensional Faddeev integral equations 
as function of the number of the grid points in the Jacobi momenta Nj ac and the spherical angles 
Ngph. The number of the grid points in polar angles is twenty. The calculations are based on the 
Bonn-B potential. 



N jac 


N sp h 


Et [MeV] 


32 


20 


-8.154 


32 


24 


-8.153 


36 


20 


-8.153 


36 


24 


-8.152 


40 


20 


-8.152 


40 


24 


-8.152 



The coupled Faddeev equations represent a set of three dimensional homogenous integral 
equations, which after discreatization turns into a huge matrix eigenvalue equation. The 
huge matrix eigenvalue equation requires an iterative solution method. We use a Lanczos- 



like scheme that is proved to be very efficient for nuclear few-body problems [33j. The 
momentum variables have to cover the interval [0, oo]. In practice we limit the intervals to 
suitable cut-offs and their values are chosen large enough to achieve cut-off independence. 
The functional behavior of the kernel of eigenvalue equation is determined by the anti- 
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symmetrized two-body t— matrix. We also solve the Lippman-Schwinger equation for the 

fully-off-shell two-body t— matrix in an approach based on a Helicity representation directly 

as a function of the Jacobi vector variables (see appendix[B]). For anti-symmetrized two-body 

t— matrix calculations forty grid points for the Jacobi momentum variables, thirty two grid 

points for the spherical angle variables and twenty grid points for the polar angle variables 

have been used respectively. Since the coupled integral equations require a very large number 

I 

of interpolations, we use the cubic Hermitian splines of Ref. [34| for its accuracy and high 
computational speed. 

TABLE VI: A list of Triton binding energy calculations ordered according to jyf x by different 
authors using slightly different numerical methods. All results for binding energies are related to 
the total isospin T = A. 



j^ ax Ref. Et [MeV] 



[35] -8.14 
36], [37] -8.17 



[38] -8.165 
[7], [22] -8.16 



[40], [41] -8.088 
[42] -8.100 



[40] -8.101 
[22] -8-103 



[43] -8.14 



32], [44] -8.13 
[7], [45] -8.14 



In Table [V] we show the convergence of the Triton binding energy as function of the 
number of the grid points for Bonn-B potential in the 3D approach. As demonstrated in 
this Table, the calculation of Triton binding energy converges to a value of E t = —8.152 
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MeV. The results of the Faddeev equations with different PW based methods are presented 
in Table [VTI in order to compare them with our calculations. The overall agreement is quite 
satisfactory. As we can see from this comparison our result provides the same accuracy while 
the numerical procedure is actually easier to implement. 



B. Expectation Value of the Hamiltonian Operator 

In this section we investigate the numerical stability of the presented algorithm and the 
3D formalism of the Faddeev equations. With the binding energy E t and the Faddeev 
component available, we are able to calculate the total wave function from Eq. ([2]) 
by considering the choice of coordinate system which is used in representation of Eq. (12 II) . 
So we can evaluate the expectation value of the Hamiltonian operator H and compare this 
value to the previously calculated binding energy of the eigenvalue equation. Explicitly we 
evaluate the following expression: 

= (tf|#o|tf> + (*|V|*> = 3 (V>|#o|*) + 3 (tf IV12I*), (34) 

where 

(V>|#o|*> = J d3p J d3q J2 J d3p ' J d3q ' (^|pq«> (pqa|#o|pVc/) (p'q'o'l*) 

-J J dx p J d(p p J dx q J d(p q 
J2^ a (pqx pq ) * a (pqx pq ), (35) 



a a' 



dpp / dqq — + — 

m Am 



x 



(*|Vi 3 |*> = J ^P f d3q Yl J d3 P J ^ (*|pqa><pqa|Vi2|p'q'a')<p'q'a'|*). 

a ' a! 

(36) 

As is well known, the rotational, parity and time-reversal invariance restricts any NN 



potential Vvi to be formed out of six independent terms 



43|, as 



Vi2(p,p') = <p|^i2| p') = ^2vi(p,p' ,x pp ,) Wi, 



(37) 



»=i 



here Vi(p scalar spin-independent functions, which depend on the magnitudes 

of the Jacobi momenta p , p' and the angle between them, x vv > = p.p', and Wi (i=l to 6) 
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are operators to the spin states of the two-nucleon such that 

Vi2m sims2 m' sim ' S2 (p,p') = (pm Sl m S2 \Vi 2 \p' m' Si m' S2 ) 

6 

= ^2vi(p,p' ,x pp/ ) {m sl m S2 \Wi\ m' Sl m' S2 ), (38) 
so the matrix elements of NN potential can be evaluated as: 

(pqalKalp'qV) = 5 3 (q - q') ( a T \V£\ a' T ) ( p a s \V 12 \ p' a' s ), (39) 

where V^ 2 is the isospin part of the potential, it is unity for the isospin-independent terms 
and T\.Ti for the isospin-dependent terms. So it can be easily evaluated as 

„ 1, isospin-independent terms: 

(a T \V l T 2 \a' T ) = T5 aTa/T , T= \ ' (40) 

I 2tf 2 — 3, isospin-dependent terms. 
The spin-space part of the potential can be evaluated as: 

(pa s \V 12 \p a' s ) = ^2^2(a s \is) (l's\a's}(Pls\Vi2\p'l's) 

7S 7s 

= ^^0^0fv<P7s|Vi2|p'Ys) 

7S 7 S 

= J2J2 9 ^y ^'7' 6m °3< 3 ( P m si m s2 \Vn\ P 1 m ' Sl m 's 2 ) 

7S 7 S 

= / ] /jQa-y 9a>Y $m S3 m> S3 V\1 m Sl m S2 m' sl m' S2 (P , p')> (41) 

7S 7s 

inserting Eqs. (jlQj) and (jUJ) into Eq. ([39]) yields: 
(pqalVialp'q'o') = 5 3 (q - q') T 5 aTQ > T ^EJE^r^V Sm * 3 < 3 V i2 m 31 m S2 m' sl m' S2 (p ,P'), 

IS 7 S 

(42) 

by these considerations the expectation value of the NN potential, Eq. f[3"6T) . can be rewritten 
as: 

(ttivi 2 i*) = EE T ^EE44^ 

/*O0 + l /*27T /*00 Z 1 "'"^ /*27T 

X / dpp 2 dx p / C% / dp p 2 / dzj, / dtp' p V 12 m Sl m S2 m>m>(p,p',Xpp>) 

Jo J-i Jo Jo J-i Jo 

/ G^g 2 / da, / dip q ^ a (p q x pq ) m a \p' q x p/g ), (43) 
Jo J -i Jo 
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X 



where x pp i = p.p' = x p x p > + — x p ^Jl — x 2 pl sin(0 p — </y) and x p > q = p'.q = x p ix q + 
v/l "4V 1 _ x 'g sin 0p' - 0?)- 

TABLE VII: The expectation values of the kinetic energy {Ho}, the NN interaction (V) and the 
Hamiltonian operator (H) calculated in the 3D scheme as a function of the number of the grid 
points in the Jacobi momenta Nj ac and the spherical angles for the Triton. The number of 
the grid points in polar angles is twenty. The calculations are based on the Bonn-B potential. 
Additionally the expectation values of the Hamiltonian operator are compared with the Triton 
binding energy results from the three-dimensional Faddeev integral equations. All energies are 
given in MeV. 



Njac 


N S ph 


(Bo) 


(V) 


(B) 


E t 


32 


20 


+39.222 


-47.356 


-8.134 


-8.154 


32 


24 


+39.222 


-47.356 


-8.134 


-8.154 


36 


20 


+39.222 


-47.357 


-8.135 


-8.153 


36 


24 


+39.222 


-47.357 


-8.135 


-8.152 


40 


20 


+39.223 


-47.358 


-8.135 


-8.152 


40 


24 


+39.223 


-47.358 


-8.135 


-8.152 



The expectation values of the kinetic energy (H ), the two-body interaction (V) and the 
Hamiltonian operator (H) are listed in Table IVlTl for Bonn-B potential calculated in the 3D 
scheme as a function of the number of the grid points in the Jacobi momenta Nj ac and the 
spherical angles N sp h- In the same Table, the Triton binding energies calculated in the 3D 
scheme are also shown in order to compare with the expectation values of the Hamiltonian 
operator. One can see that the energy expectation value and the eigenvalue energies E t 
agree with good accuracy. 

IV. SUMMARY AND OUTLOOK 

In this paper we have introduced the three-dimensional Faddeev integral equations for 
the calculation of the Triton binding energy with the spin-isospin dependent potential. In 
comparison with the PW approach, as is commonly used, this direct approach has greater 
advantages. The pertinent results can be summarized as follows: 
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1) The 3D formalism leads only to a strictly finite number of coupled three-dimensional 
integral equations to be solved, whereas in the PW case after truncation one has a set of 
finite number of coupled equations with kernels containing relatively complicated geometrical 
expressions. So the 3D formalism avoids the highly involved angular momentum algebra 
occurring for the permutations and also automatically consider all the partial waves without 
any truncation on the space part. However the 3D formulation leads to a lesser number 
of coupled integral equations in three dimensions and the PW formulations leads to more 
coupled integral equations in two dimensions. 

2) Our result for the Triton binding energy with Bonn-B potential is in good agreement 
with the pervious values calculated with the standard PW approach. The stability of present 
algorithm and the 3D formalism of Faddeev components as presented in this paper have 
been achieved with the calculation of the expectation value of the Hamiltonian operator 
and we have reached to a resonable agreement between the obtained energy eigenvalue and 
expectation value of the Hamiltonian operator. The 3N bound state calculations with AV18 
potential is also potentially valuable and the numerical results with this potential will be 
reported in the future. 

3) We predict that the incorporation of three-nucleon force probably will be less cum- 
bersome in a realistic 3D approach. This is very promising and nourishes our hope that 
four-nucleon bound state formulation and calculations with realistic two and three-nucleon 
forces in a realistic 3D approach will be more easily implemented than the traditional partial 
wave based method. 

The calculations of three-nucleon bound state, with the phenomenological Tucson- 
Melbourne (TM) 2-7T exchange three nucleon potential, and the formulation of the four- 



nucleon bound state is currently underway and they will be reported before long 



48|. 
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APPENDIX A: g ia CLEBSCH-GORDAN COEFFICIENTS 



In the usual coupling scheme, for the three identical particles with spin |, in order to 
completely classify the states of definite total spin the quantum numbers 

| 75 ) = | s 1 m Sl s 2 m S2 s 3 m S3 ) = \m Sl m S2 m S3 ), (Al) 

are replaced by the set 



| a s ) = |((si s 2 )s l2 s 3 )SM s ) = \(s 12 -)SM S ). 
The 3N basis states | as ) can be obtained from free 3N basis states \ js) as: 



(i|)|+|) = |TTt) 



(i 1)1+1) 
(i 1)1 -I) 
(i 1)1-1) 



i 

1 



III) 



fl 1)1 + I\ = _L. 



(1 |)| "I) 

(o 4)i+|> 



i 

1 



1 



ITT 
Til 



TIT 
ITI 
TTI 
ITI 



(o §)§-§) 

If one considers all total spin states, i.e. S — , 



TIT 
ITI 

TTI 
III 
TIT 
TIT 




2| ITT) 
2| Til) 



(A2) 



(A3) 



|, the relevant Clebsch-Gordan 



4 and S 

coefficients g^ a are 1, ~J 2 '~~~J 2 - ^ s indicated in section UTCl the isospin states 

are similar to the spin states, but the third component of total isospins is restricted to 
Mt = 4p for 3 He and My = +r for 3 if . Thus for a fully charge dependent calculation the 
necessary isospin coefficients g^ a are — — Since in our calculations for 

the Triton binding energy we consider only the total spin-isospin states (S — T) — (| — ~), 
therefore we only use the following Clebsch-Gordan coefficients , — W -| , , — . 



APPENDIX B: ANTI-SYMMETRIZED NN T-MATRIX AND CONNECTION 
TO HELICITY REPRESENTATION 

In our formulation of the 3N bound state, we need the physical representation of NN 
t-matrix or matrix elements a (pm Sl m S2 m tl mt 2 \t(e)\p' m' Sl m' S2 m' tl m' t2 ) a . The connection of 
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these matrix elements to those in the momentum-helicity basis is given in Ref. [30(, here we 
prepare this connection according to the notation to be used in our work. First, we introduce 
the momentum-helicity basis states for the total spin s 12 and the relative momentum p of 
the two nucleons as: 

|p;psi2A), (Bl) 

where A is the eigenvalue of the helicity operator Si2.p. By introducing parity operator P 
and the two-nucleon isospin states \ti 2 m tl2 ), the anti-symmetrized two-nucleon basis states 
are given as: 

|p; PS 12 X- t 12 y a = -i=(l - ^(-) Sl2+ ' 12 ) |i 12 ) |p; ps 12 A)„ (B2) 

with the parity eigenvalues rj n = ±1 and eigenstates |p;psi 2 A) 7r = ^(1 + VttP)\Pi P s i2^)- 
Based on these basis states the NN i-matrix element is defined as: 

^ 12 (p, p';e) ee ^(p;ps 12 X;t 12 \t(e)\p';p's 12 X , ;t 12 y a . (B3) 



As shown in Ref. 30], the selection of p' parallel to the z-axis allows, together with the 
properties of the potential, that the angular dependencies of the NN t-matrix elements can 
be simplified as: 

*aT* 12 (p, P';<0 = e" lM V ^(pv-Ae) 

= e" iM V e iA '^' t^! 2 * 12 (p, p, cos 6 pp ,- e) 

ee e i(A 'V- An »') g* u (p ) j/,cosV;E), (B4) 

the direction n pp / can be determined by the spherical and polar angles $pp/ and (f> pp >, where 

cos 6 pp i = cos 6 p cos p / + sin 9 p sin 6> p / cos(0 p — (/y ) , 
sin^pp/e^pp' = — cos 9 P sin 9 p / + sin 9 P cos 6 P ' cos((j) p — <p p ') + ism 9 p sm((f) p — (j) p /), (B5) 

and the exponential factor e 4( - A ^pp' _An ) is calculated as: 

g PP — 



In the above expressions, D s ^ 2 x ((p p 9 p 0) are the Wigner D- functions and c^ A (#) 



are rotation matrices 



31] . Finally the connection of the t-matrix elements 
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„(pm Sl m S2 m tl m t . 2 \t(e)\p' m' Sl m' S2 m' tl m' t2 ) a to those in the momentum-helicity basis, namely 
^A' 2 * 12 (P,P';£), is given as: 

a (pm Sl m S2 m tl m t2 \t(e)\p m' Si m' S2 m t m t2 ) a 

= ^ mtl+mt2 x 1+ ^ 2 e " i(Ao0p ~ A, ' ,v) E a-^(-r i2+i12 ) 

Sl2^12 



*C{~-t 12 ) mi 1 m t2 )C(^t / 12 ; m' ti m' t2 ) 
xC(--si 2 ; m, 1 m 83 )C(--ai 2 ; m' Sl m' S2 ) 

x £ OW^O^fo P'; e)- (B7) 

AA' 

It should be mentioned that ^ A J 2 * 12 (p,p', cos8 pp i; e) obeys a set of coupled Lippman- 
Schwinger equations which for S = it is a single equation but for S = 1 it is a set of two 
coupled equations (Ref. (29])- So the matrix elements of the anti-symmetrized NN t-matrix, 
which explicitly appears in Eq. (122|) . is functionally the same as Eq. flB7[) and can be 
obtained as: 



ta m si m,s 2 mt 1 mt 2 



(p,7r,x pir ;e) = a (pm Sl m S2 m tl m t2 \t(e)\nm' Si m' S2 m[ i m' t2 ) a 

= ^ mtl+mt2 ,K 1+ K 2 e " l(Ao ^ A ^ #) E ( l -^(-) sl2+tl2 : 



x C(~ -ti 2 ; rn tl mt 2 )C(--t' 12 ; m' t m' t2 ] 
x C(- -si 3 ; m Sl m S2 )C(--s' 12 ; m' Si m' S2 ) 
x E CA(^)<^(^)*ir tl2 (P 5 *l e), (B8) 

AA' 

with the same variables as Eqs. (1B5I) and (1B6I) with 7T, 0%, (f)% instead of p', 9 P >, (f) p >. 
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